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Some caveats about the evolution of the N/O 
abundance and the star formation history. 
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Abstract. We carefully analyze how the abundance of Nitrogen over Oxygen evolves when 
dependent on metallicity stellar yields with a primary component of N proceeding from 
AGBs stars are used. We show the results obtained with a chemical evolution models grid, 
calculated with variable star formation efficiencies, which produce different star formation 
histories. Finally we see how the N/O abundance is related on the evolutionary history. 
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1. Introduction 

Nitrogen abundances when compared with the 
Oxygen one, N/O, may, in principle, inform 
about the time in which low and intermediate 
mass (LIM) stars formed, and, consequently, 
when the widest episode of star formation took 
place. This idea is based on the stellar mean- 
lifetimes for LIM stars, main producers of N, 
which are longer than the ones of massive stars, 
which eject O to the interstellar medium (ISM). 

The data, however, must be carefully inter- 
preted before to reach misleading conclusions 
based in this simple scheme. If a mass of gas 
form a bulk of stars simultaneously, in a Single 
Stellar Populations (SSP), it is true that N will 
appear after O in the ISM, with a certain time 
delay. However, when data proceeding from 
different galaxies or different regions within a 
galaxy, are compared, the previous scenario is 
not longer valid, since a) most of data refers 
to spiral and irregular galaxies, where the star 
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formation is a continuous process and, due to 
that, there are a mix of SSPs; and b) the ob- 
servations correspond to the final stage after an 
evolutionary path object, and galaxies follow 
different tracks. 



2. Using the Closed Box Model 

2.1. The basic scenario 

N may be primary, NP, created directly from 
the original H by the corresponding nuclear 
reactions, or secondary, NS, if there exists a 
seed of C or O in the gas with which the star 
formed. Then the production is proportional to 
the original oxygen abundance. N is mainly 
produced as secondary in both massive and 
LIMs stars. But in LIM stars a fraction of N is 
created as primary in the third dredge-up and 
Hot B ottom Burning processes iRenzini & Volil 
(I1981I) . In massive stars N the stellar rotation 
provokes the increas e of NP, specially a t the 
lowest metallicities, (Ekstrom et al. 2008). 



2 



Molla & Gavilan: Some caveats about the N/O abundance evolution 



np+ns;- 



□ MWG HII R. 
. gal. HI! R. 

. fowZ III! ft. 
A IZ005 
. DLA 
. LMC.SMC 
SDS5 



12 + log(0/H) 

Fig. 1. The nitrogen over oxygen ratio, as 
log (N I 0), as a function of the oxygen abundance, 
as 12 + log log(0/H), for a wide sample of data as 
coded: MWG stars and HII regions; Other galaxies 
HII regions, (low metallicity objects shown sepa- 
ratel y), Damped Lyman Alpha, and SDSS estimates 
from lLiang et al.l <2006h 



Followin g the well kn own Closed Box 
Model (CBM lTinslevll 19801) . the abundance of 
metals, Z, in a region may be expressed by the 
equation: Z = plnpT 1 , where p = M gas /M tot is 
the gas fraction in that region and p is the inte- 
grated stellar yield, the elements newly created 
by a stellar generations of stars (by assuming 
that all stars that produced metals are already 
died). By applying this equation for N and O, 
and assuming that both are primary: 



Z N p N 

— = — = constant 
Z Po 
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where Zo and Z^ are the Oxygen and 
Nitrogen abundances, and po and p^j are the 
corresponding integrated stellar yields. 

If N is secondary, then: 
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We plot these theoretical lines in the dia- 
gram N/O vs O/H, Fig. Q] The first case ap- 
pears as an horizontal line while the second 
one would be represented by a straight line 
with a certain sloped When both contributions 



there exist, the horizontal line begins to in- 
crease from a certain oxygen abundance. We 
also plot the observations in this Fig. Q] (see 
references in Mol la et all 120061 It is evident 
from this diagram that both components of N 
are necessary. 

Since LIM stars produce N by both pro- 
cesses, such as mo st of authors establish (see 
iGavilan et ai]|2006l thereinafter GAV06, and 
references therein), it is immediate to try to 
explain the data by using the corresponding 
yields of LIM and massive stars. In this work 
we show how N/O evolves when a contribution 
of NP proceeding from LIM stars there exists. 

Some people think that low metallicity ob- 
jects must be young, linking wrongly oxygen 
abundance- time, and, therefore assume that 
they have not sufficient time to eject N from 
LIM stars. Because of that they search for other 
way to obtain NP. 

However, this idea that links low oxygen 
abundance to short evolutionary time is a mis- 
take: The oxygen abundance is not a time scale. 
It is possible to have a high oxygen abundance 
in a very short time, or a very low one even in 
an object which exists since a long time ago if 
it evolves very slowly dLegrandl2000f) . 

2.2. By taking into account the stellar 
mean-live-timescales 

It is possible to solve the equation of the CMB 
for both components by assuming that there are 
two classes of stars, the massive ones, m > 
8M Q , which create O and secondary N, and 
the intermediate stars, 4 < m < 8M , which 
eject prim ary N, too. The resu lting expression 
is given bv lHenrvetaHd20()0l) as: 

„ PnsPc„2 , PnsPc „3 , Pnp 

£n — — — ^ — + 



Zl+^-(Z -Z T )e->{3) 



Both lines are drawn with arbitrary yields 



2 P 2 6 P 2 o " Po 

where Z, is the abundance of each element, 
C, N, or O, p; is the corresponding yield, and 
the subscript NP and NP refer to the the com- 
ponent primary and secondary. 

Following that, N begins as secondary, 
when massive stars die, what means that N/O 
evolves over an straight line with a certain 
slope. After a delay time t, in which the oxy- 
gen abundance Z<? reaches a value Z T , the NP 
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ejected by LIM stars appears in the ISM. Just 
in that moment Z# increases exponentially, due 
to the last term of the equation, until to arrive 
to the primary N level. The behavior simulates 
a phase change, that occurs exactly when the 
first stars eject NP, time defined by the most 
massive LIM stars stellar mean-lifetimes. 

If these stars have 2, 4 or 8 M , their 
mean-lifetimes are ~ 220, 100 and 40 Myr, 
respectively, such as it is obtained from 
the functions (Fig. 13 given by the usual 
used Padova, Gen e ve ste llar tracks or from 
IBazan & Mathewsl d 19901) . The evolutionary 
tracks resulting of the above Eq. [3] are repre- 
sented in Fig. [3] panel a). 

The N/O abundance ratio increases 
abruptly when the oxygen abundance is Z r , 
which takes a different value for 2, 4, or 8 
M : The smaller the stellar mass, the higher 
the O abundance Z T at which the exponential 
function appears and the N/O increases, such 
as we show in Fig. [3^). From this plot, we 
could say that, in a Single Stellar Population 
in which NP is ejected by only a value of 
stellar mass, if, for a given O abundance, N/O 
is close to the secondary line, their stars are, in 
average, younger than those ones with a higher 
value of N/O, nearer to the primary line. 
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Fig. 2. The stellar mean-live-timescales as a func- 
tion of the stellar mass, as o btained from Padova's 
group ( P), (Bressan et al. 1993), Geneve's (G) group 
(Meynet & Maeder 2002) stellar tracks, or as given 
bv lBazan & Mathewsl Jl990h (B). 




12 + log(0/H) 

Fig. 3. The evolution of N/O, as log(N/0) vs O/H, 
as 12+log(0/H) following Eq. [3] a) for stars with 
different masses ejecting NP, as labelled and b) for 
stars of 4M Q ejecting NP but where the star forma- 
tion occurs with different efficiencies. 



2.3. The star formation efficiencies 

However, this is not the complete history. In 
order to obtain the abundance Z T , following 
Eq. Q] it is necessary to know the gas frac- 
tion, a quantity that may change with time. 
Let assume that the galaxy starts with a value 
fi = 1, what means that the mass is completely 
in gas phase, and that it decreases when stars 
forrrQ. If we assume that the gas mass g de- 
pends on time as: t| = -Hg, then g - gue~' H ', 
the mass on stars is s* = go(l - e~ ,H ), and 
^ = Tt ~ Hg, where H is the efficiency to 
form stars. In that case the star formation law 
results an exponentially decreasing function of 
the time: = Toe^'^' ', being the time-scale 
to the inverse of H. If H is high, the gas will 
consume very quickly and, therefore, the gas 
fraction /i = g/{g + s) takes a small value in 
a very short time. Obviously if H is small, the 
gas fraction maintains near from 1 . 

We show the evolution of fi(t) in in Fig.|4^) 
for some values of H as labelled. The Oxygen 



2 Note that other evolutionary scenarios are pos- 
sible 
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Fig. 4. The evolution of a) the gas fraction /i and 
b) the oxygen abundance 12 + log(0/H) for differ- 
ent values of star formation efficiency as labelled. 
The solar region value is marked as a (yellow) solid 
large dot. The solar value of the oxygen abundance 
is plotted at 8.7 Gyr by assuming that Sun born 4.5 
Gyr ago. 



abundance evolves according as shown in 
Fig. Hp). 

Over-plotted on both graphs there are 3 
lines corresponding to times, in increasing or- 
der, equal to the mean-lifetimes of stars with 
8,4 and 2 M . We obtain 3 different values 
for the oxygen abundance Z T for each evo- 
lutionary line. Equivalently, different oxygen 
abundances Z T for each stellar mass are ob- 
tained, depending on H. If we assume that stars 
of 8 M produce NP, Z T is ~ 8, 7.2, 6.8 or 
even < 6 for H= 1, 0.5, 0.172 or 00005 in a 
time of 40 Myr. That is, the NP appear in the 
ISM at an abundance that may be as lower as 
12 + log (O/H) = 6 if the efficiency to form 
stars is very low, or as higher as 8, if the effi- 
ciency is high, but always in a same time scale 
as short as 40 Myr. 

We show in panel b) of Fig.[3]the evolution 
of the relative abundance log(N/0) vs the oxy- 
gen abundance 12 + log(0/H) for three values 
of star formation efficiencies by assuming that 
NP is ejected by stars of 4 M . We see a similar 
behavior to that one shown in panel a). 



3. The metallicity dependent stellar 
yields 

The previous section results are obtained by 
taking constant effective yields p^p, Pns and 
po- Actually, the stellar yields depend on the 
metallicity Z with which stars form. And there- 
fore the effective yield for a single stellar pop- 
ulation also depends on Z. In the case of LIM 
stars, this question is even more important 
since, obviously, the proportion NP/N must 
change with Z. If the star has a very low Z, the 
abundance of O, the seed for the NS, is also 
low, so most N is created as NP. 

This is clearly seen in Fig. 2 from GAV06 
where this proportion NP/N is represented as a 
function of the stellar mass for some sets of 
stellar yields from the literatur e. In panel a) 
the so l ar abundance yields from iGavilan et all 
(12005b: Ivan den Hoek & Groenewegenldl997l) : 
Marigol d200ll hereinafter BU.VG and MA, 
respectively) are represented. In panels b), c) 
and d) the same is shown for two values of 
Z for each one of these sets separately, as la- 
belled. In all cases, NP/N is higher for the low 
Z set than for the solar abundance one. 

The dependence on Z of the integrated 
yields for a simple stellar population is shown 
in Fig. 3 from GAV06. In panel a) the inte- 
grated yield of N tends to be larger for higher 
Z. In panel b) the ratio pnp/pn clearly de- 
creases with Z, such as it must do. The differ- 
ent sets, however, have different behaviour. BU 
has smaller total yields of N than VG and more 
similar to values from MA, but the dependence 
on Z is smoother that this one shown by VG 
and MA. 

By using dependent yields on the CBM 
equation of the previous section the results 
shown in Fig [5] are obtained. Results for dif- 
ferent efficiencies H, low, solar and high, are in 
panels a), b) and c) as dotted lines. The main 
effect of changing yields with Z is that evolu- 
tionary tracks elongate when the star formation 
rate is strong, as it occurs in c). Therefore the 
efficiency to form stars, which define the shape 
of the star formation history, the level of star- 
burst with which the stars form, is essential to 
obtain an elongated or flat evolutionary track 
in the plane NO-OH. 
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Fig. 5. The evolution of log(N/0) as a function of 
12+log(0/H) for different star formation rates: a) 
low, b) intermediate and c) high efficiency, respec- 
tively. Dotted lines are the CBM results while the 
solid lines are the results obtained by GAV06 for 
different regions of the MWG disk. 



4. The chemical evolution model grid 

4.1. The local universe 

Until now only CBM results have been an- 
alyzed. However modern chemical evolution 
models are usually numerical codes that in- 
clude a larger quantity of information, taking 
into account the mean-lifetimes of stars, the 
IMF, the stellar yields, a large number of el- 
ements, etc. We have over-plotted in Fig. |5J 
with solid lines, the results from GAV06, ob- 
tained by using the multiphase chemical evolu- 
tion m odel dFerriniet alj|l994t iMolla & Dia3 
I2005I) with the VG stellar yields for different 
regions of the Milky Way Galaxy (MWG), lo- 
cated at 10, 6 and 2 kpc of galactocentric. They 
are similar to those for the CBM. 

Following the same scenario than GAV06, 
we have compute models for 44 sizes or 
masses of galaxies and 10 values of efficien- 
cies in the range [0,11 for each one of th em 
(see lMolla & DfadboOStlMolla et al.ll2006l for 
details). The results for the present time are 
shown in Fig. [6] with large open dot while the 
small (grey) dots represent the SDSS data and 
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Fig. 6. Multiphase chemical evolution model re- 
sults for the present time for a grid of 440 galaxies 
as open dots. Small dots in the background are the 
SDSS data and the other symbols are extragalactic 
and Galactic HII regions data as labelled. 



other HII regions data from the literature. Our 
results even reproduce the high dispersion ob- 
served in the plane NO-OH: there are two re- 
gions where the dispersion is higher, one lo- 
cated around (7.8, -1.45) and a second one 
around (8.7,-1), just where the low mass and 
the bright massive galaxies fall. 

These models demonstrate clearly that effi- 
ciencies to form stars have an important role in 
the evolution of the tracks in the plane NO-OH, 
and are essential to reproduce the data. 

4.2. The star formation history role 

iMallerv et all (120071) have shown that the lo- 
cation of a galaxy in the plane NO-OH is re- 
lated with the specific star formation rate, the 
star formation rate by stellar mass unit, sSFR. 
They found that the highest values of N/O cor- 
respond to the smallest values of sSFR. 

Since we also reproduce this trend, as we 
may see with different code of colors in Fig. 6, 
we may explain the subjacent reasons for this 
correlation. In our models the smallest sSFR 
occur in galaxies where the star formation was 
high in the past, in the earliest times of evolu- 
tion. The gas was rapidly consumed and there- 
fore the star formation rate decreased since 
then, showing now very small values. On the 
opposite side, when the efficiency to form stars 
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Fig. 7. a) Evolutionary tracks for six different theo- 
retical galaxies in the plane No-OH, b) Star forma- 
tion histories for the same 6 models 



is very low, the star formation history increases 
with time, and, because of that, the present star 
formation rate is high. The evolutionary track 
of an object and its final point in the NO-OH 
plane depends on the star formation history 
suffered by the galaxy. In order to demonstrate 
that, we show in Fig. [7^) 6 different evolution- 
ary tracks in the plane NO-OH and the corre- 
sponding star formation rate in panel b) of the 
same figure. 

The low mass galaxies with low efficien- 
cies have not had a high star formation in the 
past, and they may considered young from the 
point of view of the averaged age for most of 
their stars, created mainly in recent times. But 
they also have old stars able to eject nitrogen, 
that was processed as primary, since in that mo- 
ment Z was very low. 

5. Conclusions 

The present day data as well as the high- 
redshift trends in the plane O/H-N/O may be 
reproduced with chemical evolution models 
using stellar yields where LIM stars create a 
certain quantity of NP 

Differences in the star formation histories 
of galaxies are essential to reproduce the data 



and the observed dispersion, the present posi- 
tion of a galaxy in the diagram N/O vs O/H 
being determined by this evolutionary history 
It is possible to have NP ejected by LIM 
stars even at a low O abundance since O abun- 
dance is not a time scale. 
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